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QCD Calculations of Decays of Heavy Flavor Hadrons 
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Precision tests of the Standard Model and searches for New Physics in the quark flavor sector 
depend on accurate theoretical calculations of decay rates and spectra for rare, flavor-changing 
processes. The theoretical status and recent developments of techniques allowing such calculations 
are reviewed. Special attention is paid to the calculation of the B Xs^ branching ratio, the 
extraction of the fe-quark mass from a fit to _B — > Xc Ivi moments, and the determination of \Vub\ 
from spectra in the inclusive decay B — > X„ I vi . From a reanalysis of different inclusive distributions 
the updated average value \Vub\ = (3.98 ± 0.15 ± 0.30) • 10"'^ is derived. Using only the theoretically 
cleanest channels, we obtain \Vub\ = (3.70 ± 0.15 ± 0.28) ■ 10"^. 



I. INTRODUCTION 

Heavy-quark physics is a corner stone in the ongoing ef- 
fort to explore the Standard Model of elementary-particle 
interactions, to determine its parameters with greatest 
achievable precision, and to search for hints of departures 
from Standard Model predictions. Thanks to an exper- 
imental program at various facilities worldwide that has 
spanned several decades, accompanied by steady progress 
in theory, there have been tremendous accomplishments 
in this field and several important discoveries have been 
made. While the Standard Model still stands as the most 
fundamental theoretical structure explaining the funda- 
mental interactions of matter in the Universe, the com- 
bined effort of energy-frontier and luminosity-frontier ex- 
periments in the coming decade is likely to shed light on 
the deep questions concerning the physics at the tera- 
scale - questions about the origin of mass, the asymme- 
try between matter and antimatter, and the nature of 
dark matter. 

One of the main challenges in heavy-quark physics is 
to disentangle the underlying, flavor-changing couplings 
of the quarks (which in the Standard Model encode all 
information about CP violation) from the less interest- 
ing but all dominating effects of the strong interactions, 
as described by Quantum Chromodynamics (QCD). So- 
phisticated theoretical tools going beyond the realm of 
perturbation theory are required to accomplish percent- 
level precision in the calculation of decay rates and kine- 
matical distributions for the most interesting processes. 
Some of the relevant tools and recent developments are 
discussed in this talk. 



II. THEORETICAL FOUNDATIONS 



Before discussing specific applications it is important 
to summarize the theoretical tools we have for calculat- 



ing decay rates and spectra in heavy-quark physics, and 
the concepts from which these tools derive. The underly- 
ing theme is the separation of short-distance from long- 
distance physics, which is natural due to the presence of 
the large mass scale mg Aqcd , which is far above the 
scale of nonperturbative hadronic physics. Factorization 
theorems state that short- and long-distance contribu- 
tions to a given observable can be separated into Wilson 
coefficient functions Ci and nonperturbative matrix ele- 
ments Mi. Generically, up to power corrections of the 
form (AqcD/m-g)" one has 



Observable ' 



(1) 



Such a factorization formula is useful, since by virtue of it 
the dependence on the high scale mg is calculable, and of- 
ten the number of matrix elements Mi is smaller than the 
number of observables that can be expressed in the form 
shown above. Typically, this reduction is accomplished 
by means of some symmetry. The nonperturbative ma- 
trix elements can then be extracted from data or calcu- 
lated using theoretical approaches such as lattice QCD or 
QCD sum rules. The factorization scale /i in ([T]) serves 
as an auxiliary separator between the domains of short- 
and long-distance physics. Observables are formally in- 
dependent of the choice of ^; however, they inherit some 
residual dependence once the Wilson coefficients Ci are 
computed at finite order in perturbation theory. The de- 
pendence gets weaker as higher orders in the perturbative 
expansion are included. The //-independence of the right- 
hand side of ([T]) implies a renormalization-group equation 
for the functions Ci, which can be solved. In the process, 
large perturbative logarithms of the form a™ In" (rngZ/i) 
can be summed to all orders in perturbation theory. 

The formal basis of the factorization formula ([T]) is 
the (euclidean) operator product expansion. Important 
examples of this type of factorization include the effective 
weak Hamiltonian [l| 

Hcs = ^Ci{mt,mw,mz,mH,M]<!P,fi)Qi{fi) , (2) 
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which is the starting point of any calculation in quark 
flavor physics. In this case the Wilson coefRcients Ci con- 
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tain short-distance physics associated with heavy Stan- 
dard Model particles such as the top-quark, the elec- 
troweak gauge bosons, or the Higgs boson. However, if 
there exists physics beyond the Standard Model involv- 
ing some new heavy particles that couple to the Stan- 
dard Model, then the Wilson coefficients will unavoid- 
ably also be sensitive to the masses and couplings of these 
new particles. An important goal of flavor physics is to 
search for these kinds of effects. Other important ap- 
plications of the factorization theorem ^ include the 
derivation of heavy-quark symmetry relations for the ex- 
clusive semileptonic decays B D^*hDi [1, [1], and the 
calculation of the total decay rate and differential distri- 
butions in the inclusive semileptonic decay B — > Xc I vi 

Bi,i. 

Processes involving energetic light partons require a 
more sophisticated form of the factorization theorem, 
which genericaUy can be expressed as 

Observable ^ Hi{Q, ^) Ji{^/QA, ^) ® Si{fi) + ... . 

i 

(3) 

Here Q Aqcd is the hard scale of the process, e.g., 
a heavy-quark mass or the center-of-mass energy. The 
hard functions Hi capture virtual effects from quantum 
fluctuations at the hard scale. The jet functions Ji de- 
scribe the properties of the emitted coUinear particles, 
whose characteristic virtuality or invariant mass scales 
as an intermediate scale ^/QA between Q and the QCD 
scale Aqcd- Here we will only be interested in inclusive 
processes, for which the jet scale is in the perturbative 
domain {y/QA ^ Aqcd)- Finally, the soft functions Si 
describe the (nonperturbative) physics associated with 
soft radiation in the process. The symbol (E) implies a 
convolution, which arises since the jet and soft functions 
share some common variables, such as some small mo- 
mentum components of order Aqcd - These convolutions 
require distribution functions (as opposed to hadronic pa- 
rameters) as hadronic input. Examples of such functions 
are parton distribution functions, light-cone distribution 
functions, and fragmentation functions. 

The theoretical basis of the factorization theorem ^ 
is a generalization of the euclidean operator product ex- 
pansion to the time-like domain. While in the past fac- 
torization theorems such as ([3]) were derived using dia- 
grammatic methods 0, S, an elegant field-theoretic 
approach is offered by the soft-coUinear effective theory 
developed in [l^, [111) [H) [ll]- The interactions of light, 
energetic particles and the associated emissions of soft 
and coUinear radiation requires a description in terms of 
nonlocal operators, in which fields live on light-like tra- 
jectories following the paths of the energetic partons. As 
a result, this is a rather complicated effective field the- 
ory. However, once its power counting rules have been 
established, it is possible to classify operators by count- 
ing powers of Aqcd/Q- 

Many important applications of the factorization the- 
orem ^ have been discussed in the literature. In jet 
physics, factorization applies, e.g., to the cross sections 



for deep- inelastic scattering [3, [1, [1, [H, [TB], Drell-Yan 
production 0, [1, [l^, and to event shapes [13- In 
the realm of heavy-quark physics, formula ^ applies 
to the inclusive processes B I vi and B Xg'y, 

in which case the soft functions Si are referred to as 
shape functions [H, [H, [13, Uli, 111- It also applies to 
exclusive semileptonic processes such as B tt I Di and 
B — > plDi JiJ,, 23] , which have been studied in great 
detail u sing lig ht-cone QCD sum rules (for recent work, 
see e.g. [2J,[23|). Also, the QCD factorization approach 
to nonleptonic processes such as B ^ Z)(*V, B — > ttK, 
and B K*^ relies on such a factorization theorem 
[H [13, m, H [m, m. some important recent ac- 
complishments in this area have been the calculation of 
the hard-scattering kernels in the QCD factorization for- 
mula for rare hadronic decays of the form B — > hih2 
(where lu are fight pseudoscalar or vector mesons) at 
0(ag) [H, [13, [Sal, the proof of factorization for the de- 
cay B K*'y 36], and the calculation of part of the 
0{al) corrections for this process [37]. 

Any systematic theoretical approach has limitations, 
an understanding of which is necessary in order to gauge 
the accuracy of calculations and, more generally, un- 
derstand what can and what cannot be calculated in a 
model-independent way. An obvious practical limitation 
of any calculation is the need to truncate the perturbative 
expansion of the perturbative coefficients in ([T]), and 
Hi, Ji in State of the art is to compute these objects 
at next-to-next-to-leading order (NNLO) in perturbation 
theory. Even if it was possible to calculate the pertur- 
bative functions to all orders (and of course we do not 
have the tools to accomplish this!), then the asymptotic 
nature of perturbative expansions in quantum field the- 
ory would imply an irreducible ambiguity in the precision 
with which these functions can be determined, an ambi- 
guity that scales like a power of Aqcd/Q- In order to do 
better, it is necessary to include power-suppressed cor- 
rections in the calculation. In some cases such as the in- 
clusive decays B —f X^ I v and B —> Xs^, a factorization 
theorem such as ([3]) can be shown to hold at any order in 
the expansion in AQCD/'^if, fssl. [39l. |40| : however, start- 
ing already at first order in Aqcd / Tnt the number of un- 
known soft functions exceeds the number of observables. 
In other cases, such as the nonleptonic decays B — > nK 
and B — > tttt, factorization breaks down beyond the lead- 
ing power in Aqcd /™& [H, [13, [11| . It follows from these 
remarks that observables very sensitive to higher-order 
perturbative effects or power corrections cannot be pre- 
dicted with high accuracy. Well-known examples include 
some direct CP asymmetries for the decays B — > ttK and 
B — !■ TTTT, and the branching fractions for some "color- 
suppressed" decays such as B — > 7r°7r°. The fact that 
calculations are rather uncertain in these cases does not 
indicate a failure of the factorization framework (which 
is rooted in quantum field theory and as such cannot be 
questioned) , but rather that in some cases the expansion 
in powers of as and Aqcd/w;, does not converge well. 

Before turning to phenomenological applications we 
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should mention an ongoing theoretical debate on the 
question of the breakdown of factorization beyond the 
leading power in heavy-to-light transitions such as the 
_B — > TT form factor or hadronic processes such as B ^ 
TTTT. It is well known that a perturbative analysis of the 
form factor is pla gued by "endpoint divergent" convolu- 
tion integrals |4ll. l42l. I43L H . [isl ] . Power corrections to 
the factorization formula for B tttt, which arise from 
processes such as "weak annihilation", give rise to sim- 
ilar endpoint divergences [13, [11] ■ It has recently been 
claimed that these divergences can be removed, and fac- 
torization be established, in a systematic way using so- 
called "0-bin" subtractions combined with rapidity fac- 
torization Unfortunately this claim has not been 
supported by detailed calculations, and in my view it 
is unlikely that the problem can be solved following the 
approach described in that paper. Nevertheless, based 
on these results it has been argued that weak annihila- 
tion amplitudes [i^l and power-suppressed penguin am- 
plitudes [ill should be real to a good approximation, i.e., 
free of large soft rescattering phases. If true, this would 
have important implications for phenomenology. Finally, 
the prediction of the QCD factorization approach that 
charm-penguin contributions to nonleptonic decay am- 
plitudes factorize [H, [13, [ll] has been challenged in [i^ , 
and a rebuttal has been given in ^5^. More theoretical 
work on this important issue would be desirable. 

The remainder of this talk will focus on three im- 
portant applications of heavy-quark theory: the calcu- 
lation of the partial inclusive B branching ratio 
at NNLO, the extraction of the 6-quark mass from mo- 
ments of inclusive decay spectra, and the determination 
of \Vub\ from inclusive B —^ X^lvi decay distributions. 
New Physics implications oi B ^ X^^ decay have been 
discussed by G. Isidori and S. Heinemeyer at this con- 
ference, while determinations of \Vub\ from exclusive de- 
cays have been covered in the talks of E. Barberio and 
C. Davis. For a lack of time, this is a very limited se- 
lection of topics from a vast range of possibilities. My 
apologies go to all those authors whose important con- 
tributions are omitted here. 



III. THE B Xs-y BRANCHING RATIO 

The decay B Xgj is the prototype of a flavor- 
changing neutral current process. In the Standard Model 
it is mediated by loops involving heavy top quarks and 
W bosons, and it is expected that yet unknown heavy 
particles, if they exist at the TeV scale, could give con- 
tributions to the decay rate not much below the Standard 
Model level. A precise study of this process, both exper- 
imentally and theoretically, can thus serve as a precision 
test of the flavor sector and place constraints on the pa- 
rameters of New Physics models. Indeed, it is by now 
broadly appreciated that the constraints on model build- 
ing from processes such as i? ^ Xg'y can be as valuable as 
constraints from electroweak precision observables, astro- 



physical observations, and other low-energy observables 
such as electric dipole moments and the anomalous mag- 
netic moment of the muon. 

About one year ago, the calculation of most NNLO 
contributions (all that are believed to be significant) to 
the B decay rate has been completed, and a first 

prediction for the branching ratio has been derived in a 
paper by Misiak at 16 coworkers [5i|. This has been a 
heroic effort, which would not have been possible with- 
out a large collaborative approach. At the same time, a 
dedicated, model-independent analysis of the effect of a 
cut on the photon energy, > Eq with Eq in the range 
between 1.8 and 2GeV as ap plie d in the experimental 
analyses, has been performed (sl . [ssL [53 |. It is conven- 
tional (yet somewhat artificial) to write the branching 
ratio as Br(ii^^ > Eq) = Br(tot) i^(i?o), where the func- 
tion F{Eo) describes the effect of the cut, while Br(tot) 
is the "total" branching ratio (which in practice refers 
to a very low cut E^ > 1 GeV, since the total branching 
ratio does not exist). 

We begin with a brief review of the NNLO calculation 
of the total branching ratio. The first task was to com- 
pute the Wilson coefficients in the effective weak Hamil- 
tonian ^ at NNLO. This required two- loop matching 
calculations for the coefficients Ci_6 at the weak scale 
[55] and three-loop matching calculations for the coeffi- 
cients and Csg [56\ . These matching conditions are 
the place where New Ph ysic s contributions would enter 
the theoretical analysis CT]. Next, the Wilson coeffi- 
cients have to be evolved down to a scale /Lt ~ us- 
ing the renormalization group. As input, one needs the 
three-loop anomalous dimension matrices for the opera- 
tors Qi-e [III and Q7J, Qsg [El], as well as the four-loop 
anomalous dimensions accounting for the mixing of Qi-e 
into Qj-y, Qsg [60,]- With these results at hand, the second 
task was to compute the matrix elements of the various 
operators at NNLO. The full two- loop corrections to the 
matrix element of the leading dipole operator Q-^j were 
obtained in [bJ, [13, [H, [i^l ■ For the required three- loop 
matrix elements of the operators Q1.2 only an extrapola- 
tion in the mass of the charm-quark is available based on 
calculations performed in the artificial limit rric 3> »n&/2 
[65| . One of the dominant uncertainties in the prediction 
results from this extrapolation. The remaining matrix 
elements are known at 0(/3oq;^) [H, [l^l- The resulting 
prediction for the total branching ratio in the Standard 
Model is [U 

Br(£;^ > 1 GeV) = (3.27 ± 0.23) • 10"^ . (4) 

The combined theoretical uncertainty results from a 3% 
error due to higher-order perturbative corrections, a 3% 
error in the extrapolation in the charm-quark mass, an- 
other 3% uncertainty from parameter variations (from 
TOt, rric, Us etc.), and finally a 5% error assigned to ac- 
count for nonperturbative power corrections to the total 
rate. 

The result ([5]) is a mile stone in the history of heavy- 
flavor theory. Note that the largest contribution to the 
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theoretical uncertainty is now due to our lack of control 
over nonperturbative power corrections. Unfortunately, 
it will be very difficult if not impossible to reduce this er- 
ror with current technology. A recent analysis has identi- 
fied a new class of nonlocal power corrections to the total 
decay rate, which are described in terms of matrix ele- 
ments of trilocal light-cone operators ^6§\. Naive model 
estimates of these effects point to a small reduction of the 
total rate with an uncertainty of 5%. However, this num- 
ber is only a guess, and it is difficult to substantiate it. 
There are other, similar types of nonlocal contributions, 
which are currently being analyzed [69j . 

The presence of a cut > Eq on the photon energy, 
which is applied in order to reduce the background in all 
measurements of the B — > Xsj branching ratio to date, 
leads to significant complications in the theoretical anal- 
ysis. The reason is that it introduces a sensitivity to the 
low scale A = mi, — 2Eo « 1 GeV, which is barely in the 
perturbative domain. One thus needs to deal with a com- 
plicated many-scale problem characterized by the hierar- 
chy nib > aA^^A > a > Aqcd • A systematic framework 
for this situation is the multi-scale operator product ex- 
pansion developed in [tO, [zH ■ In this case the theoretical 
accuracy that can be achieved is not set by size of the ex- 
pansion parameters as(mb) « 0.2 and AQCD/»7ib ~ 0.1, 
but rather by as{A) « 0.4 and Aqcd/A « 0.5 (though 
first-order corrections in this ratio are absent). A detailed 
analysis at NNLO comes to the conclusion that [H^l 

F(1.6GeV) = 0.93i°:°^, (5) 

which tends to be a stronger suppression than that pre- 
dicted by naive model estimates [72[ . The above number 
includes a perturbative uncertainty of %, parameter 
dependences of 2%, and an error from power corrections 
of also 2%. 

Combining the results Q and ([5]), one obtains the 
NNLO prediction ^ 

Bt{E^ > 1.6 GeV) ^ (2.98 ± 0.26) • 10"'' . (6) 

It is consistent within errors with the result (3.15±0.23) • 
10"'' quoted in [5l|, which is based on a naive pertur- 
bative expansion in agimb). For completeness we also 
quote the number (3.49 ± 0.49) • lO"'' obtained in ^ 
using a renormalon-inspired model for the shape func- 
tion in B ^ Xs"f decay. Note that the difference with 
respect to ([B]) is primarily due to a larger total branching 
ratio obtained in this scheme, not so much due to the 
treatment of cut-related effects. 

The theoretical result ^ for the branching ratio is 
about 1.4 standard deviations lower than the current 
world-average experimental value Br(_E-y > 1.6 GeV) = 
(3.55±0.26)-10-* [zl], indicating that there IS some room 
for possible New Physics effects. We emphasize, however, 
that to date no actual measurement has been performed 
with a cut as low as 1.6 GeV. Indeed, the quoted exper- 
imental number involves a model-dependent extrapola- 
tion to low energy. In the future, it would be desirable 



if experiment and theory would be confronted at a value 
of Eq that can indeed be achieved experimentally. The 
lowest such value at present is Eq — 1.8 GeV. 



IV. EXTRACTION OF mt FROM A FIT TO 
Xclvi MOMENTS 

The operator product expansion allows for a model- 
independent calculation of moments of decay distribu- 
tions in the inclusive semileptonic process B — s- Xc I 

Q . Hadronic physics is encoded in a few parameters 
defined in terms of forward i?-meson matrix elements of 
local operators. They are called /i^ (or Ai), ^Xq (or A2), 
etc. The most important parameters in the prediction 
are, however, the heavy-quark masses vrib and TOc. The 
only assumption underlying the theoretical calculations 
is that of quark-hadron duality. It is believed to be reli- 
able for the relevant energy release Ai? ~ — m^i. 

Traditionally, the parameters { | Vcb \ , mb , mc , /i J , /i^ } 
are extracted from a global fit to experimental data on 
moments of the lepton energy and invariant hadronic 
mass spectra in B ^ Xc I vi decay and of the photon 
energy spectrum in i? ^ Xs^ decay. The experimental 
data include measurements reported by the BaBar, Belle, 
CLEO, CDF, and DELPHI collaborations [74], while the 
theoretical calculations are based on formulae derived in 
[ziSHIlil- The global fit takes account of the strong 
correlations between the various quantities. The status of 
the theoretical calculations is such that the leading terms 
in the operator product expansion for the moments are 
know at 0{(}Qa1) but not yet at 0{a1). Power-suppressed 
contributions are only known at tree level. It will be 
important in the future to increase the accuracy of the 
calculations. The technology for obtaining the complete 
0(0;^) corrections to the leading term exists [t^, and it 
is expected that results will soon be published. Also, 
the 0{as) corrections to the leading power-suppressed 
effects, those proportional to the parameters /zj and /z^., 
are currently being calculated. Results have already been 
published for the terms proportional to /i^ [s^l , but they 
have not yet been included in the fits. 

The fit strategy is as follows: A set of experimental re- 
sults, including errors and correlations, is fitted to a set of 
theoretical equations derived using the operator product 
expansion. The fit is performed in the space of param- 
eters {jVcbl, 7716, TOc, /J.^, mentioned above. Besides 
\Vcb\, these parameters must be defined in a particular 
renormalization and subtraction scheme. The most pop- 
ular schemes are the kinetic scheme 18 ill, the IS scheme 
[S^ . and the shape- function scheme |23. Issj . While in 
principle it is merely a matter of choice which of these 
schemes one prefers, fitting in scheme A and convert- 
ing the results to scheme B might give slightly differ- 
ent results from fitting directly in scheme B, if trun- 
cated perturbative expressions are employed. This fact 
should be taken into account when deriving estimates 
for the theoretical uncertainties; specifically, the quoted 
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TABLE I: Translation of two different mj, values from the kinetic scheme to other schemes, at different orders in perturbation 
theory. The errors in columns two to four indicate the scale dependence (m(,/2 < /i < 2m(,). The last column shows some 
reference values for comparison. 



m^'" = 4613 MeV 


0{as) 


0(/3oa?) 


o(«') 


Reference values 


nibirnb) [MeV] 
mf [MeV] 
[MeV] 


4317 
4698 1 1 
4609 1 2 


4159 
47451 I 
4584 


41951^1 
47421^1 

4643 lit 


4164 ± 25 (e+e- hadrons [84]) 
4701 ± 30 (fit in IS scheme [74]) 
4630 ± 60 (used by HFAG [74]) 


m^'" = 4677 MeV 


mt{mb) [MeV] 
mf [MeV] 
mf^ [MeV] 


437711^ 
4761 1 I 
46711 2 


4217 1;^^ 

48071 9 
4648 


42531^^ 
48041^1 
4706 111 


4164 ± 25 {e+e- ^ hadrons [84]) 
4751 ± 58 (fit in IS scheme [74]) 
4707 ± 56 (HFAG update for LP07) 



uncertainties should not be less than the uncertainties 
inherent in the scheme translations. To illustrate this 
point, we quote results derived by the Heavy Flavor Av- 
eraging Group (HFAG) for this conference. They ob- 
tain m^'" = (4.613 ± 0.035) GeV from a fit performed 
in the kinetic scheme and = (4.701 ± 0.030) GeV 
from a fit performed in the IS scheme [zl]. On the other 
hand, translating the first number from the kinetic to 
the IS scheme at OiPoal) would give m^^ from kin ^ 
(4.745 ± 0.036) GeV, which is about 1.5 standard devi- 
ations higher than the value obtained from the fit in the 
IS scheme. This indicates that the uncertainty is larger 
than the 30 MeV quoted from the fit. 

Table [J collects results obtained by translating the cen- 
tral values m)j'" = 4.613 GeV (upper portion) and m^'" = 
4.677 GeV (lower portion) from the kinetic scheme into 
other schemes (MS, IS, and shape-function scheme). The 
formulae translating from one scheme to another involve 
perturbative expansions, and following common practice 
we choose to evaluate the coupling as in these expansions 
at a scale fi varied between mf,/2 and 2mi,. This results 
in the quoted uncertainties. A few comments are worth 
mentioning: First, note that the scale uncertainty from 
the conversion into the MS scheme remains significant 
even at two-loop order. This unfortunate fact prevents 
us from using the very precise value for mf,(mf,) recently 
obtained from an analysis of the total e+e" hadrons 
cross section in the threshold region [s^j to derive simi- 
larly precise values for the 6-quark mass in other schemes. 
Next, note that the value of nib at NLO in the IS scheme 
has a tiny scale uncertainty, yet the value obtained at 
NNLO differs from the NLO value by more than 40 MeV. 
This illustrates that scale variation can sometimes under- 
estimate the true theoretical uncertainty. Finally, note 
that for the translation into the shape-function scheme 
two-loop corrections not associated with running cou- 
pling effects appear to be important. It would therefore 
be desirable to include the full 0{aj.) corrections in the 
moment fit. 

A serious problem with the way the fit is currently per- 



formed by HFAG is that besides the moments of semilep- 
tonic spectra in _B ^ Xc I vi decay also moments of the 
B — !■ Xgj photon energy spectrum are used. Based on 
what we know about B — > Xgj decay, this treatment can 
no longer be justified theoretically! Whereas a clean the- 
oretical approach to the B —^ Xc I vi moments is provided 
by the operator product expansion, the theoretical basis 
for the calculation of the B Xg'^ moments is far more 
complicated and subject to uncontrollable theoretical un- 
certainties. First, it is unquestionable that the existing 
measurements of the B Xsj moments are performed 
in a kinematic region where shape-function effects are 
still important [85|. To correct for these effects, a (shape- 
function) model-dependent "bias correction" is applied 
to the data [s^. More generally, it has now been estab- 
lished that as soon as operators other than the dipole 
operator Q-^j are not ignored, then there is no local oper- 
ator product expansion for the B — > Xgj decay rate and 
moments, even outside the shape- function region [68l[69j. 
Indeed, a novel factorization formula can be derived for 
the photon energy spectrum, which is graphically rep- 
resented in Figure [TJ Besides a "direct photon" contri- 
bution (first term), there appear "single resolved" and 
"double resolved photon" contributions, which are im- 
possible to calculate with present tools. Their impact on 
the moments of the photon spectrum is currently under 
investigation. 

When the goal is to extract heavy-quark parameters 
such as rub in a model-independent way, one should thus 
not include the B Xgj moments in the fit. An analysis 
based only on the B ^ Xclvi moments leads to m\™ = 
(4.677 ± 0.053) GeV in the kinetic scheme and = 
(4.751 ± 0.058) GeV in the IS scheme tZ4]. The increase 
of the uncertainties can be understood by observing that 
in the combined fit there was s slight tension between the 
results from the B — > Xc I vi and B Xg^ moments, as 
illustrated in Figure [2l The increased errors shown above 
more properly reflect the true uncertainties. These errors 
will be reduced (perhaps to 30-40 MeV) once the 0{a1) 
and 0{as/m^) corrections will be included in the fit. 
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FIG. 1: Graphical illustration of the QCD factorization theorem for B Xs'y decay in the endpoint region [63|- The dashed 
lines represent soft interactions, which must be factored off the remaining building blocks to prove factorization. 
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0.2 I — ' — ' — ' — ' — ' — ' — ' ' — ' — ' — ' — ' — ' — I 

4.4 4.5 4.6 4.7 

m, (GeV) 

FIG. 2: Results for the parameters mi, and (defined in the 
kinetic scheme) obtained from moment fits to B — > Xc I vi 
(blue), B Xs'y (green), and both (red). The plot is an 
update of a figure in J2j (courtesy of HFAG). 



When the new values for rub are translated into the 
shape-function scheme, one obtains m^^ = (4.707 ± 
0.056) GeV, which is larger than the value (4.63 ± 
0.06) GeV used so far by HFAG. This observation has 
important implications for the determination of \Vub\- 

V. INCLUSIVE DETERMINATION OF Kf,| 

The theoretical basis for the analysis of differential 
decay distributions in the inclusive semileptonic process 
B Xu iDi is a. factorization formula of the form ([3]), in 
which the soft functions Si are generalized parton distri- 
bution functions for the B meson called shape functions 
[11, [MHO, HH HI. The leading-order shape function is a 
universal, process-independent quantity, which also (to a 



large extent) determines the shape of the B Xsj pho- 
ton energy spectrum. The strategy is therefore to extract 
the shape function from the B ~^ Xs^ photon spectrum 
and then use this information to predict B I 
decay distributions [2^, ^7, ^] . The functional form of 
the shape function is constrained by moment relations, 
which relate weighted integrals over the shape function to 
the heavy-quark parameters nib, A ^tc. extracted from 
the B Xc I yi moment fit. An alternative is to employ 
shape-function independent relations between weighted 
B Xsi and B XJvi spectra [l^, HI, S IM HI • 
Both approaches are equivalent; yet, not all predictions 
that have been obtained using the shape-function inde- 
pendent relations are up to the standard of present-day 
calculations. 

The elimination of background events from B — > 
Xcli'i, in which the charm quark is misidcntified, is 
accomplished by means of different experimental cuts. 
The most common ones are a cut Mx < to_d on the 
hadronic invariant mass of the final state, a cut Ei > 
(m^ — rnjj)/2mB on the energy of the charged lepton, 
a cut > (ms — mu)'^ on the invariant mass squared 
of the lepton- neutrino pair, or a cut < mjj/niB on 
the plus component of the total momentum of the final- 
state hadrons. All are designed such that they eliminate 
charm background while keeping some portion of the sig- 
nal events. From a theoretical perspective the ideal cut is 
that on hadronic invariant mass, followed by the P_|_ cut. 
They keep more than 50% of the signal events. Cutting 
on lepton energy or leptonic mass is far less efficient. 

Until this summer, the most complete theoretical anal- 
ysis of inclusive B X„ / Di spectra was based on calcu- 
lations described in and referred to as the BLNP 
approach. It includes complete perturbative calculations 
at NLO with Sudakov resummation, subleading shape 
functions at tree level, and kinematical power correc- 
tions at 0{as). The calculation of 0{al) corrections is in 
progress and could have a significant impact [1^ . An al- 
ternative scheme called "Dressed Gluon Exponentiation" 
(DGE) [93, [H, [9^ employs a renormalon-inspired model 
for the leading shape function, which is less flexible in its 
functional form that the forms used by BLNP. Also, no 
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TABLE II: Compilation of \Vub\ values obtained using different partial B — ^ Xuli>i decay rates, analyzed using the BLNP 
approach [23, H^. The first column of values is obtained by extracting the fe-quark mass from the combined fit to 

B Xc I vi and B — > Xs^ moments [tJ ■ The second column is obtained by only using the theoretically clean information 
from B — » Xc I vi moments (courtesy of HFAG). 



Metnod 


a 


\Vub\ [10 J (old J 


\Vub\ [10 \ (new) 


mf [GeV]: 




4.63 ±0.06 


4.71 ±0.06 


CLEO, El > 2.1 GeV 


UA 


3.91 ± 0.46 ± 0.44 


3.52 ±0.41 ±0.35 


Belle, El > 1.9 GeV 


9.5 


4.67 ±0.43 ±0.36 


4.35 ±0.40 ±0.33 


BaBar, Ei > 2.0 GeV 


11.4 


4.23 ±0.24 ±0.39 


3.89 ±0.22 ±0.33 


BaBar, Ei > 2.0 GeV, s^="' < 3.5 GeV^ 


13.8 


4.37 ±0.29 ±0.49 


3.94 ±0.27 ±0.39 


Belle, Mx < 1.7 GeV 


10.5 


3.92 ±0.26 ±0.32 


3.66 ±0.24 ±0.27 


BaBar, Mx < 1.55 GeV 


13.6 


4.09 ±0.20 ±0.39 


3.74 ±0.18 ±0.31 


Belle, Mx < 1.7 GeV, > SGeV^ 


9.3 


4.23 ±0.45 ±0.36 


3.97 ±0.42 ±0.31 


Average 




4.31 ±0.17 ±0.35 


3.98 ±0.15 ±0.30 



attempt is made to include subleading shape functions, 
which among other things means that the predictions 
for more inclusive partial rates would not be in accord 
with the operator product expansion beyond the leading 
power in Aqcb /n^b- The DGE model is therefore less 
rigorous than the BLNP approach, even though it leads 
to numerical results that are compatible with those of 
BLNP. Very recently^ Gambino et al. (GGOU) ^ have 
produced a code for B — > I vi spectra that is compara- 
ble in its completeness with that of BLNP. Working in the 
kinetic scheme, these authors use moment constraints to 
model the subleading shape functions differently for the 
various structures in the hadronic tensor. They include 
the recently calculated 0{(3Qa'j.) corrections to the decay 
rates [l^ , which were not yet available for BLNP. On the 
other hand, no attempt is made to resum Sudakov loga- 
rithms. The numerical results obtained using the GGOU 
code are consistent with those of BLNP in both their 
central values and error estimates. 

In Table HIl we compile results for \Vub\ obtained from 
various various experimental cuts. The data have been 
analyzed using the BLNP approach. It is crucial to re- 
alize that these results depend very sensitively on the 
shape function used in the analysis, and in particular on 
its first moment, which is related to the 6-quark mass. 
The coefficient a given for each cut indicates how the 
partial rate scales with rrifc. The extracted value of \Vub\ 
thus scales like (mb)~°/^. The third column in the table 
shows results reported previously by HFAG, which are 
based on using the old value m^^ ~ (4.63 ±0.06) GeV for 
the 6-quark mass. Following my request, HFAG has now 
redone the analysis using the new value m^^ — (4.71 ± 
0.06) GeV obtained without including the B Xg'y mo- 
ments in the global fit. The new combined average value 
is 



CLEO (EJ 






3.52 ±0.41 ±0.35 ' ' " 






BELLE sim. ann. (m,^, q") 


A 




3.97 ±0.42 ±0.31 






BELLE (EJ 






4.35 ± 0.40 ± 0.33 






BABAR (EJ 






3.89 ±0.22 ±0.33 






BABAR (E^, s;;~) 






3.94 ±0.27 ±0.39 






BELLE (mx) 






3.66 ±0.24 ±0.27 






BABAR (m^) 






3.74 ±0.18 ±0.31 ' ' 






Average +/- exp +/- (mb,theory) 






3.98 ±0.15 ±0.30 






x'/tiof=6.3/6(CL= 39%) 
OPE-HQET^SCET (BLNP) 






Phys.Rev.D72:073006,2005 
mi, input from b— > c 1 V moments 

III 




HFAG 




LP 2007 
1 1 1 


2 


4 





iKbl = (3.98 ± 0.15 ± 0.30) • 10" 



(7) 



FIG. 3: Compilation of \Vub\ determinations from inclusive 
spectra in B — > Xu I vi decay, using = (4.71 ± 0.06) GeV 
for the 6-quark mass, as determined using only model- 
independent information for B Xc I i^i moments (courtesy 
of HFAG). 



It is significantly lower than the previous value. A graph- 
ical representation of the individual determinations is 
shown in Figure [3] Note that there is rather good agree- 
ment among the different measurements. From a theo- 
retical perspective, one would give preference to the ex- 
tractions based on the most efficient cuts, such as the 
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cuts on Mx or P+ (not reported in the table). For these 
quark-hadron duahty is expected to work best and, more 
importantly, the potential contamination from weak an- 
nihilation effects is minimal. Averaging the BaBar and 
Belle results derived form the Mx cut, my personal "best 
guess" for |Kb| is 

\Vub\ = (3.70 ± 0.15 ± 0.28) • 10"^ . (8) 

The most important conclusion from this discussion is 
this: When only model-independent information is used 
to determine the 6-quark mass and thereby derive con- 
straints on the leading shape function, then the value 
extracted for \Vub\ from inclusive decays is in rather 
good agreement with that derived from exclusive de- 
cays. For instance, recent analyses of the distribu- 
tion in__B TT Ivi decays using unquenched lattice QCD 
[99I llOOj ] or light-cone QCD sum rules [11] yield values 
in the range between (3.3-3.6) • 10^'', with uncertainties 
of about lg'4 • 10~^ f74|. In my view there is no obsta- 
cle to averaging the "inclusive" and "exclusive" results 
to obtain a single value for \Vub\, but this is best done by 
HFAG. 



VI. CONCLUSIONS 

An important theoretical challenge in heavy-flavor 
physics is to disentangle strong-interaction effects from 
the underlying weak-interaction transitions, which probe 
the flavor structure of the Standard Model and mea- 



sure fundamental parameters such as the elements of the 
CKM matrix. Tremendous progress in theory has been 
made based on heavy-quark expansions and factorization 
techniques combined with effective field-theory methods. 
Thanks to these developments, and to the hard work of 
many theorists, we now have accurate predictions, in- 
cluding high-order perturbative effects and nonperturba- 
tive power corrections, for a number of important decay 
processes, including in particular the inclusive processes 
B Xsj and B Xc^ulvi. The "|T4b| vs. sin 2/3 puz- 
zle" , i.e., the fact that there used to be a tension between 
the determination of \ Vub\ from inclusive semilcptonic de- 
cays as compared with the indirect value inferred from 
the measurement of the angle [3 of the unitarity trian- 
gle, is resolved when only model-independent informa- 
tion is used in the extraction of the 6-quark mass and 
other heavy-quark parameters. The new value of \Vub\ 
obtained in this way is in good agreement with the de- 
termination from exclusive decays. 
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